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1. Introduction 

Binary alloys in stoichiometric compositions invariably exhibit ordered structures at 
low temperatures. As we depart from perfect stoichiometric compositions, it is not 
possible to populate the lattice in the given compositions so as to produce a perfectly 
ordered structure. Take for example, a B75A25 binary alloy on an fee lattice. One 
of the possible stable ordered phases is the L12 arrangement as shown in figure 1. 
An example of this is CU3AU. In a cubic unit cell the corner is occupied by an A 
atom, while the three face centres by B atoms. Since there are N corners and 3N face 
centres (N being the total number of unit cells in the solid) and exactly as many A 
and B atoms, at this composition the L12 ordered arrangement exactly populates all 
the lattice sites. Ordered arrangement becomes impossible, say for B 7 oA 30 . However, 
the following arrangement becomes possible : since there are now 1.2N A atoms, N 
of them may occupy the N corners. The 3N face centres may be occupied randomly 
by the 2.8N B atoms and 0.2N remaining A atoms. The original A sublattice remains 
ordered, while the B sublattice becomes disordered. Since there are on the whole 3N 
face centres, the occupation probability of the A and B atoms in this sublattice are 0.93 
and 0.06 respectively. This arrangement is quite different from the completely random 
alloy, where all sites are randomly occupied by either the A or B atom with probabilities 
0.3 and 0.7 respectively. It is also rather different from the partial disorder defined by 
T0|1 . In this communication we shall define partial disorder in the manner described 



above. 



The recursion method was introduced by || |7j as a convenient and numerically efficient 
method for calculating Green functions and realistic physical properties like the local 
density of states (LDOS), the Fermi energy and the band energy. The method comes 
to its own in situations where translation symmetry of the potential in an effective one- 




B sublattice 



A sublattice 



Figure 1. The L12 atomic arrangement for 75-25 AB binary alloys on a fee lattice 
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electron Hamiltonian is lost. The Bloch Theorem is violated, and we cannot employ the 
standard reciprocal space techniques unless we carry out a homogeneity restoring mean- 
field approximation (like the coherent potential approximation, CPA). This can happen, 
for instance, at a surface or in a random alloy. In particular, if the ion-cores sit on a 
topologically distorted network or at a rough surface or interface, where the recursion 
method seems the only natural way of dealing with electronic property calculations. 

There exists a large body of literature on the application of the recursion method 
to various situations ([§). For substitutional disorder, the recursion method has 
been employed in tandem with the augmented space formalism (|L1|, [12], [13], [14], [T5|] ) 



(augmented space recursion or ASR) for dealing with random alloys, taking into account 
local configuration fluctuations like chemical clustering, local lattice distortions due to 
size effects and inhomogeneous randomness at a surface because of surface segregation. 
In spite of this extensive body of literature, the recursion method does not seem to have 
gained wide acceptability in the electronic structure community, and one often has had 
to face skepticism regarding its accuracy and feasibility when applied with sophisticated 
electronic structure techniques to realistic systems. 

In this communication we shall first present the ASR in detail for the most general case 
of systems with many atoms per unit cell. Our subsequent aim will be to study partially 
disordered systems in alloys with off-stoichiometric compositions. 

We shall carry out recursion calculations based on the tight-binding linearized muffin-tin 
orbitals (TB-LMTO) method ([1], 0, [§) for various metals and alloys, to start with for 
perfectly ordered systems, compare them with standard reciprocal space methods also 
based on the TB-LMTO and then carry out a thorough analysis of the convergence and 
accuracy of the technique. The aim will be to have a convincing handle on the possible 
errors in the method with a view to minimizing them within our tolerance limits. Having 
convinced ourselves about the feasibility of the method, we shall first apply it to a toy 
model for partial disorder before carrying on to realistic systems. 

2. Formalism 

2.1. Generalized Augmented Space Recursion within the TB-LMTO formalism. 



In earlier works ( [JT3J, [14|, [L5|]) it has been established that for a disordered system the 
Augmented Space Theorem maps a disordered Hamiltonian described in a Hilbert space 
H onto an ordered Hamiltonian in an enlarged space \l/, where the space \l/ is constructed 
by augmenting the configuration space $ of the random variables of the disordered 
Hamiltonian together with the Hilbert space H of the disordered Hamiltonian . Thus 
^ = $ ® H. The configuration average of the Green function reduces to the evaluation 
of a particular matrix element of the resolvent of this enlarged Hamiltonian in the 
augmented space. Hence, if one performs a recursion in the augmented space, one can 
obtain the matrix element necessary to calculate the configuration average of a Green 
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function directly. The advantage of this method is that it does not involve a single site 
approximation or the solution of any self-consistent equation (which is a prerequisite 
for the CPA or its generalizations). Furthermore, one can treat both diagonal and off- 
diagonal disorder on an equal footing. In spite of its immense potential, the method 
could not be implemented in a real alloy systems successfully earlier because of the 
large rank of the augmented Hamiltonian. For a binary alloy with N sites and with only 
s-orbitals it is Nx2 N ). However, if we exploit the symmetry of the augmented space 
(which arises due to the homogeneity of the disorder) the rank of the irreducible part 
of the augmented space in which the recursion is effectively carried out, is considerably 
reduced and recursion becomes tractable. Unlike other methods like the embedded- 
cluster method, where the configuration averaging is done by explicit averaging over 
all the distinct configurations of the embedded cluster (144 for a fee lattice with first 
shell of neighbours) , this generalized augmented-space recursion technique yields, in a 
single recursion, the configuration average directly. In a recent paper, || has rederived 
the augmented space technique and related it to the idea of spatial ergodicity which 
connects volume averaging to configuration averaging in the infinite solid limit. 

As has been mentioned earlier, since the recursion method needs a localized, short- 
ranged basis for its operation, one can implement augmented space recursion in 
the framework of the TB-LMTO formalism. We now describe the methodology of 
generalized augmented space recursion in the framework of the TB-LMTO. 

The second order TB-LMTO Hamiltonian is written in terms of potential parameters 
and screened structure matrix. For a random binary alloy A^B^, the LMTO Hamiltonian 
in the most localized representation is given by, 

H {2) = E u + h-hoh (1) 

where, 



h — 2_^ [C RLa — E RLa J VrLo + Z^ Z^ ^RLa SRLa,R'L'a' ^-R'L'a 1 ^RLa,R'L'a' (2) 

RLa RLa R'L'a' 

C, o and A are potential parameters of the TB-LMTO method, these are diagonal 
matrices in the angular momentum indices. Also o~ l has the dimension of energy and is 
a measure of the energy window around E around which the approximate Hamiltonian 
7i 2 is reliable. 
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R denotes a cell position label associated with a TB-LMTO basis and L = (£mm s ) is the 
composite angular momentum index, a denotes an atom in the i?-th cell whose position 
is R + C, a - tIr is the site-occupation variable which takes values or 1 depending upon 
whether site a in the i?-th cell is occupied by an A or a B atom. For partial disorder 
this is a random variable whose probability density depends upon which sublattice it 
belongs to, hence the label a associated with it. The structure matrix §RLa,R'L'a' is 
non-random in case of substitutional alloys with negligible size mismatch. Now one can 
obtain the full 7i 2 by inserting h in expression (1). 

Vrlol and Trl<x,r' v a! are the projection and transfer operators in Hilbert space H 
spanned by tight-binding basis {\RLa}}. 



Vrlu = \RLa)(RLa\ 

TRLa,R>L>a> = \RLa) (R'L'a'\ (3) 

The expanded Hamiltonian 7i in the augmented space is constructed by replacing the 
random site occupation variable n^ by its corresponding operator representation Ai^ 
in configuration space. 



h = E ((Cf " E ul) Z + ($C L - SE uL ) M%) ® V RLa + . . 



RLa 



RLa R'L'a 



where, Ai% is given by 



M a R = x a A VL + x a B VL + y/(x a A x%){r^ + T^ 

Now inserting equation(4) in equation(l) one can obtain the full Hamiltonian 7i This 
Hamiltonian is now an operator in a much enlarged space \1/ = $ <g> 7i (the augmented 
space). The Hilbert space \& is spanned by the basis set {\R,L,a,{C})}, where {C} 
is a pattern (TTIT • • •}) spanning the configuration space and described in terms of the 
cardinality sequence {C} of positions of the j in the pattern. The enlarged Hamiltonian 
does not involve any random variables but incorporates within itself the full statistical 
information about the random occupation variables. 
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In the presence of off-diagonal disorder, which is invariably present in the form of the 
TB-LMTO Hamiltonian, even the reduction of the rank of the invariant subspace on 
which recursion acts, using the symmetries of the augmented space, does not allow us 
to sample as many configuration states as we would like. This is because, as recursion 
proceeds, the number of configuration states sampled at the n-th step of recursion 
becomes unmanageably large when n > 5 (for example). To do away with this problem, 
the working equations are transformed so as to put the Hamiltonian for the recursion 
in an effective diagonal disorder form. This allows one to sample further shells in 
augmented-space and to confirm the shell-convergence of the recursion. 

To do this, we first suppress all the indeces and write the expression for the resolvent 
as follows : 
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Using the augmented space theorem, we can write the expression of configuration 
averaged Green function as, 



<GW*La(£)> =(R,L,a,m\(EI-H) 1 |i?,L,a,{0}) 
where, 

A a L (A" 1 / 2 ) \R, L, a <g> {0}) + Fl (A" 1 / 2 ) \R, L,a® {R}) = |1} 

1/9 

We define [A" (1/A)] ' |1} as the normalized initial recursion vector |1). 
Finally, we arrive at the convenient expression : 

< G RLa , RLa (E) >= (1| [E - A + B + F - S + (J + S) o (j + S)]' 1 11) 
where, 
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here, J = Ja + Jb + Jf and 6 = 04 + 6b + dp, where, 
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where, 



A2(Z)=x3Zf + xSZf 
Bt(Z)= (x%-x a A ) (Zi-Zl 



FZ{Z) = ^x\x% [Zt-Zi 
Z is any single site parameter. 

Though the computational burden is considerably reduced due to diagonal formulation, 
the recursion now becomes energy dependent as is clear from the expressions of B and JF 
above. The old formalism was free of this constraint. This energy dependence makes the 
recursion technique rather unsuitable because now we have to carry out one recursion 
per energy point of interest. This problem is tackled using seed recursion technique 
([|J). The idea is to choose a few seed points across the energy spectrum uniformly, 
carry out recursion on those points and then fit the coefficients of recursion throughout 
the whole spectrum. This way one can save huge computation time. For example, if 
one is interested in an energy spectrum of 100 points, in the bare diagonal formulation 
recursion has to be carried out at all the 100 points but in the seed recursion technique 
one needs to perform recursions only at 15-20 points. The whole idea stems from the 
fact that in most of the cases of binary alloys, it is seen that the recursion coefficients 
a n and (3 n vary quite weakly across the energy spectrum. So, one can easily pick up a 
few of them and fit throughout the whole range of energy by a suitable function. 

2.2. Convergence of the recursion method 

Before we apply the above methodology to study realistic alloy systems, we must first 
ensure that the technique is stable and convergent. We shall illustrate the nature of 
convergence of various calculated physical quantities which are basic for any electronic 
structure calculation. This study of convergence is essential because otherwise we will 
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Figure 2. The Density of States for Cu and CU3A11 using (top) the recursion method 
and (bottom) the k-space method 



not be able to set the cut-offs of various parameters like number of shells in real and 
configuration spaces, recursion levels and so on, for further calculations and ensure that 
our numerical results are within our tolerance window. 

Before we set out to analyze the errors in the method, let us first compare the results for 
the density of states for Cu, CU3AU, V and Ni as obtained from a forty-step recursion 
carried out on an exact twenty shell real-space map. The choice of the systems is 
deliberate. Cu is a full rf-shell noble metal, V is a half-filled metal in the lower end 
of the transition series while Ni is a magnetic transition metal. CU3AU is an ordered 
alloy which is stable in the L12 arrangement on a fee lattice. Figures |2| and |^ show this 
comparison. 



For Cu and Cu 3 Au most of the detailed features of the density of states calculated from 
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Figure 3. The Density of States for Ni and V using (top) the recursion method and 
(bottom) the k-space method 



the recursion calculations match with the k-space results. As expected, the k-space 
results show sharper structures. The Fermi energies in both cases match to within a 
few hundredths of a rydberg. Both the majority and minority band density of states 
for ferromagnetic Ni are reproduced excellently well in the recursion calculations. The 
agreement is equally good for V upto the Fermi energy. The unoccupied part of the band 
shows disagreement. However, it must be understood that the recursion is done with the 
second order Hamiltonian in the most localized representation of the TB-LMTO, while 
the k-space calculations are done with the orthogonal representation. A higher order 
recursion calculation should improve the unoccupied part of the band considerably. 

How good are our results ? When we talk of convergence of the recursion method, 
we have to be careful in stating precisely what we mean. Finite space approximants 
to Green functions do not converge for real energy values. This problem arises in 
every computational method, as noted by || . The problem definitely arises in the 
usual k-space integration techniques, where methods using either complex energies or 
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complex k-s have been attempted. The cause of this non-convergence is that an arbitrary 
small perturbation, like adding a single atom to a large but finite system, can shift all 
eigenvalues of the system. This causes an infinite change in the Green function near 
its corresponding poles. Thus, the precise meaning of the convergence of the recursion 
should imply rather the convergence of physical quantities built out of it. Most physical 
quantities are averages over the spectrum of the type : 



pE F 

F{E) = / f{E')n{E')dE' 

JEn 



Eq is the lower band edge, and f(E) is any smooth, well-behaved function of E. It is the 
convergence of these quantities which will decide whether the recursion is convergent or 
not. For example, the Fermi energy is defined by 



/ F n(E')dE' = n £ 

J-Ea 



where n e is the total number of electrons. While the band energy is 

/E F 
E' n(E')dE' 
-E 

We shall study, in general, the convergence of indefinite integrals of the kind 
M k (E) = f {E') k n{E')dE' 

J — OO 

fk 

The integrand E is monotonic and well behaved within the integration range. A 
measure of the root-mean-square error in the moments is 

A < = {e^e„ (C WW™ - (S E E ? mE)dE )) } (5) 

Errors can arise in the recursion procedure because of two distinct sources : (i) the error 
that arises because we carry out a finite number of recursion steps and then terminate 
the continued fraction using one of the available terminators ; (ii) the error that arises 
because we choose a large but finite part of the nearest-neighbour map and ignore the 
part of the augmented space very far from the starting state. Haydock has justified both 
these approximations by stating that (i) if the terminator is so chosen so as to reflect the 
asymptotic behaviour of the continued fraction, errors should be small, and (ii) since 
we can write expressions for the continued fraction in terms of self-avoiding walks in 
augmented space, long walks are dominated by those that wind round the starting state 
and do not go far away from it. 

We shall first carry out a simple error analysis of the continued fraction expression for 
the Green function because of errors created on the continued fraction coefficients. 
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The recursion is a two-term recurrence relation. We may therefore generate from this, 
a pair of linearly independent set of polynomials through the relations : 

bn+\X n+ i{E) — (E — a n )X n {E) — 6„X n _i 

where, X n (E) is either P n (E) or Q n (E) according to the initial conditions : 

P 1 (E) = 1 P 2 (E) = (E - ai)/&2 
Qi{E) = Q 2 (E) = 1 

The approximated Green function in terms of the terminator T(E) is given by : 

Q N+1 (E) - b N Q N (E)T(E) 



G(E) 



h [P N+1 (E) - b N P N T(E)} 



The terminator determines entirely the essential singularities of the the spectrum, [g] 
showed that a finite composition of fractional linear transformations like the one above 
can at most add a finite number of poles to the spectrum. The essential singularities of 
the exact G(E) and T(E) coincide. The fractional linear transformation redistributes 
the spectral weights over the spectrum. 

Let us now assume that we make errors {Sa n , Sb n } in the corresponding continued 
fraction coefficients. If we now start generating the orthogonal polynomials, starting 
from the exact initial conditions, but with the errors in the continued fraction 
coefficients, we shall obtain a pair of sets {P n } and {Q n } ■ In general we shall have, 

P n (E) = (1 + A n (E)) P n (E) + B n (E)Q n (E) 

If we substitute this back into the recurrence relation and keep only the first order terms 
in the errors, 



A n (E) = {Sa n [P n+1 (E)Q n+1 (E)} 

+5b n [P n (E)Q n+1 (E) + P n+1 (E)Q n (E)]} /h 
B n (E) = {-5a n P n+1 (Ef - 6b n [2P n (E)P n+1 (E)}}/( 



Using the above and the expression for the local density of states , we find that the first 
order relative error produced in the local density of states 
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Figure 4. The root mean square errors in the first (full lines), second (dashed lines) 
and third (dashed dotted lines) moments of the density of states as functions of the 
recursion steps 



where R(E) = TZe G(E). If we define the weighted Hilbert transforms of P n (E) as the 
so-called associated functions : 



Q n (E) = ne\[_ 



P n {E')n{E') 
E-E' 



dE' 



These associated functions are also solutions of the three-term recursion. They are not 
polynomials, but are nevertheless orthogonal to the set P n (E). 

In terms of these, the error in the density of states is : 
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(6) 
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If the first N continued fraction coefficients are exact, that is in case we carry out our 
recursion on a N-shell neighbour map upto N steps and then terminate, the error in the 
various moment functions are : 



oo 

SM k (E) = +J2{6a n A£\E) + 6b n +iB$°\E)} 



n=N 



Where, 
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P n+1 (E')Q n+1 (E')(E') k n(E')dE' 
P n+1 (E')Q n+1 (E')(E') k n(E')dE' 



From this expression and equation(6) we can obtain an expression for the overall error 
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Figure 6. The convergence with recursion step of the band energies of Cu,Cu3Au, Ni 
and V 



in the moments. Numerical results for the errors in the moment functions are shown in 
figure [|. The convergence in Fermi and band energies with number of recursion steps 
are shown in figures || and ||. 

These are calculated on a 20 shell neighbour map and upto 40 steps in recursion. The 
figures clearly illustrate the convergence of the procedure with increasing recursion steps, 
as surmised by Haydock. As expected the higher moments converge more rapidly. The 
rapid convergence of higher moments is at the basis of the reproduction of most of the 
density of states shape in an approximate recursion procedure. It should be noted, 
however, that we carry out recursion beyond 40 steps, very soon the procedure becomes 
unstable. This instability arises due to two sources : (i) we have carried out recursion 
on a 20 shell neighbour-map. As we go beyond 20 steps of recursions, finite size effects 
begin to show up. Numerically till 40 steps these errors are tolerable. Beyond 45 they 
lead to instability. We can control this by increasing the size of the neighbour-map. (ii) 
Numerical cumulative errors lead to the loss of orthogonality in the recursion generated 
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Figure 7. The convergence with recursion step of the magnetic moment of Ni 

basis. This can be controlled by deliberately orthogonalizing after every 20 steps or so. 
In the presence of disorder, the sharp structures in the density of states are smoothened 
by finite life-time effects. The convergence of the higher moments is even more rapid. 
The above analysis in ordered systems is therefore in a worst scenario situation. 

Figure 7 shows the convergence of the magnetic moment of Ni as a function of the 
number of recursion steps towards the value 0.6 [iB obtained from fc-space techniques. 
Magnetic moment is a physical quantity which is very sensitive to the errors in numerical 
approximations. The figure shows that within 35-40 recursion steps on a 20 shell real- 
space map, the magnetic moment has converged within our error tolerance. 

It is evident from the above analysis that the size of the neighbour map and the number 
of recursion steps required can vary from system to system. For every situation, we 
have to carry out this error analysis before we can rely on our numerical results with a 
degree of satisfaction. 

3. Calculations on a model system 



Before we carry out calculations on a real alloy, let us first apply our formalism to a 
toy model in order to understand the effects of partial disordering. We shall consider a 
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Figure 8. Ordered atomic arrangement on a square lattice for a A 50 B 5 o alloys 

50-50 AB alloy ordered first on a square lattice as shown in figure |8| and only s-states 
on this atomic arrangement. The Hamiltonian in a tight-binding basis set is then, 
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where {ij} denotes that i and j are nearest neighbour cells on the lattice, a is 1 or 2 
according to whether we a referring to the corner atoms or central atoms in a square 
unit cell (dark and light atoms in figure [§D. If we consider only the nearest neighbour 
overlaps at a distance a/\r2 where a is the square lattice constant, the diagonal and 
off-diagonal terms of the Hamiltonian are (referring to figure ^j 



€ia,ja' 



t 



ia,ja' 



t 



ia,ja' 




Ti-Tj are in the (10) and (01) directions 



Ti-Tj are in the (10) and (01) directions 



The ordered lattice has a density of states which has a central band gap with integrable 
infinite Van Hove singularities at the two internal band edges, two kink singularities 
within the band and the usual square-root singularities at the external band edges. As 
expected, the density of states is symmetric about the band centre e = 0.5. This is 
shown in figure |TI]. 

The density of states for the perfectly disordered lattice, where each lattice site is 
occupied by either A or B atom with probabilities proportional to their concentrations, 



is shown in the figure |TT|. Disorder washes away the Van Hove singularities and the 
central band gap is filled up. The results are identical to a disordered square lattice 
with lattice vectors (a/y/2,0) and (0,a/y/2). 
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t=1 nearest neighbour at ~ 0.7 a 

t«1 next nearest neighbour at a 

Figure 9. Nearest neighbour overlaps for a central cell in a square lattice with lattice 
constant a units 




Figure 10. Density of states for the perfectly ordered lattice 



The next figure [12] shows the density of states and Fermi energies for the partially 
disordered half-filled alloys at just off-stoichiometric compositions (49.5-50.5) and (50.5- 
49.5). Disorder washes out the Van Hove singularities, although vestiges of the kink 
singularities remain. The signature of the internal infinite singularities show up as peaks, 
but disorder fills up the internal band gap. The band-edge square-root singularities 
remain as artifacts of the termination procedure. Loss of stoichiometry weights the two 
'bands' differently leading to a loss of symmetry about the band centre. 
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Figure 11. Partial Density of states (states/energy-cell) for the A and B 
atoms (top, bottom) for a completely disordered alloy at different off-stoichiometric 
compositions 



The figure [13| shows the band energy for the partially ordered alloys. There is a jump 
at the stoichiometric composition with the band energy for the ordered alloy falling half 
way between the two branches. This jump can be understood by carefully examining the 



graphs on the left panels and right panels of the figure |13[ The band energy represents 
the average energy (centre of gravity) of the portion of the graph to the left of the 
Fermi-energy The density of states on the left panel have greater weight below the 
Fermi energy as compared to those on the right panel, ensuring the band energy to be 
more negative for the compositions sown on the left panels. The sudden jump is the 
result of the Fermi energy shifting across the pseudo band-gap as we cross the 50-50 
composition and including the high peak at E=0, which shifts the average energy to 
higher values. The specific behaviour of the band energy depends sensitively on the 
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Figure 12. Density of states (states/energy-cell) for partially disordered alloy, at 
different off-stoichiometric compositions 



features of the density of states. The specific behaviour for our toy model may not 
obtain for realistic alloy systems. 

In contrast, the behaviour of the band energy with composition for the fully disordered 
alloy is smooth, reflecting the smooth behaviour of the density of states without any 
internal band gaps. The figure indicates that for compositions to the left of the 50-50 
stoichiometric one, the alloy prefers to be partially ordered. At 50-50 the ordered alloy is 
stable for our toy model. While for compositions to the right of 50-50 the alloy stabilizes 
in the completely disordered phase. 

4. Remarks and Conclusion 



In this communication we have first proposed a generalization of the ASR for partially 
ordered binary alloy systems. Partial ordering of the type where disorder in different 
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Figure 13. Band energies for the partially ordered, fully ordered and fully disordered 
alloys at stoichiometric (50-50) and nearby off-stoichiometric compositions 



sublattices are different has been studied in particular. We have analyzed the 
convergence of the recursion technique in the worst scenario case of no disorder and 
have argued that disorder smoothens out structure in the density of states, so that for 
disordered alloys the convergence of the moments is even faster. Finally we have studied 
a model system, in order to understand the effects of partial disordering. We are now 
in a position to apply our formalism to realistic alloy systems. This will be the subject 
matter for a subsequent communication. 
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